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We discuss the ability of the SKA to observe QSO proper motions induced by long-wavelength gravitational 
radiation. We find that the SKA, configured for VLBI with multiple beams at high frequency (8 GHz), is sensitive 
to a dimensionless characteristic strain of roughly 10"^"', comparable to (and with very different errors than) other 
methods in the 1/yr frequency band such as pulsar timing. 



1. Introduction 

The history of astronomy has been that of our 
measurement of photons from astrophysical ob- 
jects and, more recently, the early universe (the 
CMB). Analogously, the universe is thought to 
be suffused by gravitational radiation — a back- 
ground from sources in the early universe (e.g., in- 
flation) and "foregrounds" from the gravitational 
interaction of massive astrophysical objects, such 
as the massive black holes in the centers of galax- 
ies. This gravitational radiation perturbs space- 
time, altering path-lengths and lensing light rays, 
throughout the universe. The small amplitude of 
the perturbations (strains of /i ~ 10~^^) makes 
these effects very difficult to observe. In this sec- 
tion, we examine the time-dependent "gravita- 
tional lensing" which results in apparent proper 
motions of distant galaxies, potentially observ- 
able by the SKA. 



2. GW-induced Proper motions 

The pattern of proper motions induced by a 
single gravity wave with "-I-" polarization, am- 
plitude h, and frequency uj in the z direction is 



- Luh sin cut s'm 9 {6 cos 2(1) ~ (/)sin2(/)) . 

(1) 



The pattern from an arbitrary wave can be cal- 
culated by superposition of this pattern at dif- 
ferent amplitudes and directions. Obviously, the 
particular pattern is not preserved by superpo- 
sition, but the power spectrum of the pattern is 
Pj. We calculate the power spectrum by decom- 
posing the pattern into vector spherical har- 
monics, Y^(i), {A = E,B) the generalization 
of the usual scalar functions Y^™ to vectors de- 
fined on the sphere. In fact, these are closely 
analogous to the tensor spherical harmonics used 
in the study of weak lensing and of the polariza- 
tion of the CMB. Indeed the goal - to measure 
the power spectrum of the geometrical pattern - 
is much the same in both cases. We decompose 
the total pattern into 
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where the coefficients af^ are functions of fre- 
quency and the aim separates out the spatial and 
frequency dependence. The quantity hc{w) is the 
characteristic strain at a given frequency We 
can finally define the power spectra 
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However, we do not directly observe the proper 
motion, /i, itself, which is an angular velocity. 
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Rather, we measure the position of objects at two 
(or more) epochs separated by time St - the inte- 
gral of the angular velocity. The power spectrum 
of this position difference is then an integral over 
the power spectrum. In analogy to Eq. 12 we 
define 9im,{St) as the transform of the position 
difference pattern, with power spectrum 



= 2 



du; 



M/(cj) [1 - cos{ujSt)] .(4) 



In fact, Gwinn et al |SI show that a full 5/6 
of the power in gravitational radiation is at £ = 
2 (the quadrupole); in practice there is little 
need to measure the entire spectrum. However, 
the iime-frequency (w) spectrum depends on the 
frequency content of the gravitational radiation 
background: 



and 
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where r2Gw(^) gives the contribution to the cos- 
mological critical density of the GW background 
per logarithmic integral of frequency. 

3. Proper motion error analysis 

From the study of the CMB, we recall that the 
error at a particular multipole is given by |5l()j 
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where w gives the weight (inverse variance) per 
solid angle and we have assumed observations uni- 
formly spread around a fraction /sky of the sky. 
This formula is just the usual one for the variance 
of the square of a quantity with known variance: 
{2£ + l)/sky is the number of modes at a given £. 
The Ml term then gives the sample variance, and 
the term the noise variance (i.e., resulting 
from interferometric phase errors). 

It remains, then, only to insert the expected 
signal for Mi and the noise characteristics of the 
SKA. In the band of frequencies probed by the 
SKA, the dominant contribution is likely to be 



massive black holes (MBHs) in close binaries at 
the centers of galaxies: 
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where M^^^ = AfiM2(Mi +M2)"^/^ gives the so- 
called "chirp mass" of the system, angle brackets 
denote the average over the massive black hole 
mass function, and the dimensionless factor Ifi2 
gives an integral over the merger rate of black 
holes lij. 

Beneath this foreground potentially lies the 
cosmological background from an epoch of infla- 
tion jZ|: 
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where is the energy scale of inflation, nx is 
the Tensor spectral index from inflation (usually, 
riT — and nx < 0), and we have ignored the 
so-called running of the index; this formula gives 
a generous upper bound and thus any realistic 
background is well below any astrophysical fore- 
grounds such as Eq. |S| 

From Eqs. we see that we need integrals 
of the form 



dujhl{uj) [1 - cos{ujSt)] . 



(10) 



For these power-law spectra, we of course just 
pick out the amplitude ai lu ^ 1/St. 

4. Discussion 

For the proposed measurements with the SKA, 
we will likely only be placing limits on the grav- 
itational radiation background and so the noise 
term will be dominant. The weight per solid an- 
gle is given hy w — Na^"^ /{^irfs^^y) for N proper 
motion observations with error tr^ on each ob- 
servation. From Ref. we can expect an rms 
position accuracy of lO/ias = 5 x lO^^^rad with 



the SKA used in VLBI mode at (say) 8 GHz. 8. Fomalont & Reid, this volume. 

Without "muhibeams," the SKA will be able to 9. Jarvis, M., Rawlings, S., this volume, 

measure the positions of roughly 10^ point-like 
QSOs in about one month; with 100 multi-beams 
this increases to the roughly 10^ objects avail- 
able on the sky [5]. This gives a total weight 
= (10^^^)^ for the multibeam configuration 
and /sky = 1/2. This indicates that our sensitiv- 
ity to gravitational radiation will be at a level of 
roughly h ~ 10~^^, competitive with other meth- 
ods such as pulsar timing 21- In rnore detail, we 
see that this noise contribution will likely still be 
dominant for ui ~ yr~^ for the signals in Eqs.lHUHI 
Without multi-beaming, we of course lose another 
factor of 10 in sensitivity to h, and the method 
is less competitive. However, this calculation as- 
sumes that the position errors are uncorrelated; 
in practice phase errors from the troposphere will 
correlate position errors, and this can mimic the 
low angular frequency pattern of the gravitational 
waves. We have also ignored the loss of some low- 
frequency information to the "geoedetic" infor- 
mation required for VLBI, as discussed in Ref. [3]. 
Nonetheless, these are very different systematic 
problems than those to which pulsar timing will 
be sensitive. 

We thus propose to use the SKA configured for 
VLBI in multibeam mode at 8 GHz, observing a 
preselected population of 10^ pointlike AGN in 
two epochs, separated by at least one year (more 
epochs and a longer time difference both increase 
the sensitivity). 
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